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ABSTRACT 

The  ocean  and  atmosphere  arc  driven  by  the  fluxes  of  momentum,  heat,  and 
water  mass.  The  importance  of  the  fluxes  of  momentum  and  heal  is  well  recog¬ 
nized  by  both  meteorologists  and  oceanographers.  The  hydrological  cycle  is, 
however,  only  given  considerable  attention  in  atmospheric  models  since  the  latent 
heat  release  is  an  important  source  or  the  atmospheric  general  circulation.  The 
hydrological  cycle  is  given  less  attention  in  ocean  models  although  it  is  realized 
that  evaporation  and  precipitation  arc  contributors  to  the  surface  buoyancy  flux 
which  determines  the  depth  of  mixing  and  drives  the  thcrmohnlinc  circulation. 

The  cloud  and  the  ocean  mixed  layer  arc  highly  coupled  by  both  the  heal  and 
moisture  fluxes  across  the  air-ocean  interface.  Two  lime  scales  arc  found  in  this 
paper:  a  sea  surface  temperature  (SST)  lime  scale,  r7  ,  that  is  virtually  controlled 
by  the  oceanic  planetary  boundary  layer  (OPBL),  and  a  cloud-SST  coupling  time 
scale,  t„  r  .  These  two  time  scales  depend  on  the  stability  of  the  marine  atmo¬ 
spheric  boundary  layer  (MABL). 

An  air-ocean  model  is  developed  in  this  paper  for  the  coupled  system  of  an 
unstable  atmosphere  overlying  a  stable  ocean.  The  model  results  demonstrate 
that  the  exchanges  of  heal  and  water  fluxes  across  the  sea  surface  leads  to  both 
growing  and  decaying  modes  of  oscillation  on  3-6  day  and  20-30  day  lime  scales. 
These  oscillatory  solutions  arc  entirely  thermodynamic  and  do  not  require  wave 
dynamics  for  their  existence. 

1.  Introduction 

Since  1970  significant  progress  has  been  made  both  in  our  ability  to  carry  out 
air-sea  interaction  field  work  and  in  our  understanding  of  many  of  the  processes 
found  on  both  sides  of  the  air-sea  interface.  Many  studies,  both  observational  and 
theoretical,  have  shown  that  the  surface  wind  and  the  SST  arc  two  important  el¬ 
ements  in  the  air-sea  coupled  system  (Chu,  1989).  Since  clouds  have  significant 
effects  on  the  large-scale  atmospheric  circulation  through  the  transfer  of  heat  and 
moisture,  and  on  the  OPBL  through  the  attenuation  of  the  solar  radiation  at  the 
ocean  surface,  there  may  be  strong  feedback  between  the  clouds  and  the  OPBL. 
Fig. I  shows  the  main  physical  processes  (heat,  mass,  and  momentum  fluxes)  at 
the  two  adjacent  boundary  layers:  the  OPBL  and  the  marine  atmospheric 
boundary  layer  (MARL).  Fig. 2  illustrates  the  feedback  pathways  between  the 
clouds  and  the  ocean  mixed  layer. 
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Fig. I  Main  physical  processes  al  the  two  adjacent  boundary  layers. 


Fig. 2  Feedback  paths  between  OPBL  and  clouds. 


The  feedback  mechanism  between  clouds  and  ocean  mixed  layer  in  the  buoy¬ 
ant  forcing  regime  (upward  buoyancy  flux  at  the  ocean  surface)  was  illustrated 
by  Chu  and  Garwood  ( 1 989).  However,  the  TOGA  (Tropical  Ocean  Global  At¬ 
mosphere  Program)  observations  have  shown  the  warm  pool  regions  of  the  west¬ 
ern  Pacific  to  be  in  the  buoyant  damping  regime  (downward  buoyancy  flux)  for 
the  OPBL  due  to  fresh  water  influx  caused  by  an  excess  of  precipitation  over 
evaporation.  The  OPBL  is  usually  stable,  and  the  MABL  is  generally  unstable. 
Therefore,  the  feedback  mechanism  in  this  coupled  system  (unstable  MABL  - 
stable  OPBL)  should  have  different  form  from  the  buoyant  forcing  ease. 

Although  our  coupled  model  is  one-dimensional,  we  arc  aware  of  the  impor¬ 
tance  of  the  horizontal  advcction  and  the  limitation  of  one-dimensional  models. 
However,  the  intent  of  this  work  is  to  develop  a  formalism  to  examine 
thermodynamic  feedback  between  the  two  fluids.  Because  we  wish  to  concentrate 
on  the  thermodynamic  interaction,  horizontal  advcction  is  ignored  initially. 
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2.  Similarity  Functions  in  the  MABL 

According  to  the  similarity  theory  it  is  assumed  that  if  a  variable  (wind,  tem¬ 
perature,  or  moisture,  etc.)  is  appropriately  sealed,  then  its  profile  follows  a  uni¬ 
versal  function  whose  form,  in  general,  is  determined  empirically.  For  the 
barolropic  atmosphere  the  appropriate  scales  for  wind  and  height  arc 


u  =  K 


au„ 


I /I 


(I) 


where  Ux  is  the  gcostrophic  wind  speed,  which  is  assumed  here  to  he  I ()/;//.?  ,  u„. 
is  the  atmospheric  friction  velocity.  The  gcostrophic  drag  coefficient  C  is  defined 
by  Yamadn  (1976)  as 
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(2) 


where  k  is  the  von  Karman  constant,  0.4.  The  roughness  parameter,  is  ap¬ 
proximately  I.SxIOVn.  The  ratio  hjr^  is  0.6  xIO7  for  an  MABL  height. 
Iis  =  I  km.  The  heat  and  moisture  transfer  coefficients  arc  determined  hy 
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where  w,  0,  q  arc  the  vertical  velocity,  potential  temperature,  and  specific  humid¬ 
ity,  respectively.  The  superscript  '  '  “  means  the  disturbance,  and  '*  '  indicates 
the  values  being  taken  at  the  top  of  the  MABL.  The  heat  transfer  coefficient  is 
further  computed  by  Yamnda  (1976)  as 

=  =  (4) 

where  Pr0  is  the  turbulent  Prandtl  number  for  neutral  stability,  having  a  value  of 
0.74  according  to  Busingcr  cl  al.  (1971).  Based  on  the  Wangara  data,  the  simi¬ 
larity  functions  A,  B,  C  arc  experimentally  determined  by  Yamnda  (1976)  for  the 
stable  atmosphere: 

A  =  I  .855  -  0.38 hJLa  (0  £  /»„//.„  <;  35), 

A  =  -2.94 (HJLm  -  I9.94),/2  (35  <  hJLm)\ 

B  =  3.02  +  0.3 hjl.a  (0  £  hJLa  <  35). 

B  -  2.85(/itf/La  -  1 2.47)1'2  (35  <  HJLJ; 

C  =  3.665  -0.829 (0  £  hJLa  <,  1 8), 

C=  -4.23  (hJLa-  II.2I)1'2  (IX</ra/ZJ  (5n) 

and  for  the  unstable  atmosphere  (hJLm<  0) 

A  -  10.0  -S.  145(1  -  O.OOX376/i„/ 
n  =  3.02(1  - 


C  =  1 2.0  -  8.335(  I  -  0.03 1 0 OhJLJ  1  /3  (5/») 

where  /.„  Is  the  atmospheric  Obukhov  length  scale.  The  computation  for  the 
moisture  transfer  coefficient  is  not  very  clear  yet.  In  this  paper  we  assume  that 

Q  =  Q/  (b) 

Substitution  of  (5  a,b)  into  (2)  and  (4)  leads  to  the  apparently  strong  dependence 
of  Cv  C„.  f|  on  the  atmospheric  stability  parameter  hJL„,  as  shown  in  Fig. 3. 
These  parameters  have  much  larger  values  for  the  unstable  than  for  the  stable 
atmosphere,  i.c., 

Cg  ~  0.03 1 6,  C„.  CF  ~  0.063  {hjl.a  <  0) 

Cf  —  0.003,  C,h  CF  ~  0.004  (hJLa  >  0)  (7) 


K 

Fig.3  Impendence  of  Cg,  C,„  C,  on  the  atmospheric  stability 
parameter  ( hjl *). 

X  Time  Rale  of  Change  of  the  Cloud  Cover 

The  time  rate  of  change  of  cloud  cover  is  assumed  to  lie  proportional  to  the 
moisture  supply  divided  by  the  amount  of  water  vapor  necessary  to  produce  the 
model  cloud.  The  main  processes  causing  the  cloud  dissolution  arc  precipitation 
and  mixing  with  the  environmental  air.  The  cloud  evaporation  due  to  mixing  with 
ambient  air  is  a  complicated  problem,  and  is  neglected  for  the  sake  of  simplicity 
here.  Thus  the  equation  for  cloud  cover  is 

i)n  _  (*',o  T  £  —  r,) 

iit  hf  {) 

where  ht  is  the  total  amount  of  water  vapor  needed  to  create  the  cloud  over  a  unit 
area.  From  mean  distributions  of  temperature  and  mixing  ratio  in  the  environ¬ 
mental  air  outside  the  cloud  and  inside  a  deep  cumulus  cloud  (Kuo,  1963),  we 
estimate  that  cm.  The  large-scale  horizontal  moisture  convergence  in  the 

column  of  atmosphere  per  unit  area  is  denoted 
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4.  Relationship  between  Precipitation  Rate  and  Cloud  Cover 

By  linear  regression  of  hourly  rain  amounts  and  satellite  IR  brightness  data 
obtained  during  Phases  I,  II,  and  III  of  GATE,  Albright  ct  al.  (1985)  presented 
a  linear  relationship  between  average  precipitation  rate  P,  in  boxes  1 .5*  (I  68  km) 
on  a  side  and  cloud  cover  n  of  the  boxes  by  clouds  with  lops  colder  than  -36"  C: 

=  (0.472  +  8.333  n)  x  |0-7  (9) 

This  result  confirms  Arkin's  (1979)  earlier  analysis  for  the  GATE  B-scalc  array. 


5.  Ocean  Mixed  Layer 


Much  of  the  one-dimensional  theory  for  the  OPBL  or  mixed  layer  is  depend¬ 
ent  upon  the  validity  of  two  crucial  hypotheses.  The  first  of  these  is  that  vertical 
mixing  within  the  turbulent  boundary  layer  and  cnlraimcnt  mixing  at  its  base 
occur  in  response  to  the  local  atmospheric  forcing  -  the  surface  wind  stress  and 
the  buoyancy  flux  at  the  sea  surface.  The  second  hypothesis  is  that  the  mechan¬ 
ical  energy  budget  is  the  key  to  the  understanding  and  prediction  of  mixed  layer 
dynamics  (Garwood,  1977,  1979). 

The  buoyancy  flux  is  attributable  l<>  heal  flux,  evaporation,  and  precipitation. 
The  shear  production  of  turbulence  is  attributable  to  surface  wind  stress.  In  the 
western  Pacific  warm  pool  regions,  the  excess  precipitation  effect  prevails  over  the 
buoyancy  effect  of  heat  lost  at  the  ocean  surface,  causing  the  net  buoyant  flux  to 
be  downward.  The  mixed  layer  depth  equals  the  Obukhov  length  scale: 
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2*V 
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(10) 


where  L*  is  the  oceanic  Obukhov  length  scale.  The  oceanic  friction  velocity  is 


/  ^  ,1/2 
■  I  p  ) 


and  B  is  the  downward  surface  buoyancy  flux 


/?  =  /!*.?( />,-£)- 


(II) 


where  ps,p9  arc  the  air  and  sea  water  densities,  a  is  the  sea  water  thermal  ex¬ 
pansion  coefficient,  and  ft  is  the  salinity  contraction  coefficient.  Because  the 
OPBL  is  in  the  buoyant  damping  regime,  the  mixed  layer  heat  and  salinity 
budgets  include  no  entrainment: 


fit 
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(12) 


0S_  =  (E  W 

dt  liw 


(13) 


where  is  the  specific  heal  under  constant  pressure.  The  parameters  A,  and  As 
.arc  the  horizontal  advcclion  for  temperature  and  salinity,  respectively. 


The  cUccls  of  clouds  on  the  buoyancy  flux  at  the  ocean  suifacc  arc  two-fold: 
(I)  decreasing  through  the  increase  in  the  net  heal  loss  at  the  ocean  surface,  F 
,  by  reducing  the  incoming  solar  radiation,  and  (2)  increasing  due  to  precipi¬ 
tation.  The  net  heat  loss  from  the  ocean  surface  is  given  by 

r«Rh-R,+  -  0 ]  -I  LpJCJJ'fifaSJJ  -  $]  (14) 

where  q,( T)  is  the  saturated  mixing  ratio.  The  incoming  solar  radiation  absorbed 
by  the  ocean  surface  is  R„  and  Ru  is  the  net  energy  loss  from  the  ocean  surface 
by  longwave  radiation,  L  is  the  latent  heat  of  vaporization  of  water,  and  //,  is  the 
sensible  heat  flux  to  the  air.  The  standard  bulk  formulae  arc  used  to  calculate  the 
surface  evaporation: 

e = pjcjUfAiATj  -  q]lnw  ( 1 3) 


6.  Cloud  Effects  on  the  Net  Radiation  at  the  Ocean  Surface 

Clouds  reduce  the  solar  radiation  incident  at  the  ocean  surface  by  scattering 
and  absorption,  which  is  computed  by  Budyko's  (1978)  formula 

=  [I  -  ~  «.,«( I  -  '»)]K,n  ( I <*) 

Here  Rrt)  (340  W  m  J)  is  the  solar  radiation  absorbed  by  the  ocean  surface  layer 
under  a  clear  sky.  The  parameters  a,„  and  a,n  represent  albedos  of  the  earth- 
atmosphere  system  with  complete  cloud  cover  and  a  cloudless  sky,  respectively, 
and  have  the  following  values: 

asn  ~  0.46,  otin  =  0.2 

The  ocean  surface  emits  longwave  radiation  to  the  atmosphere  and  to  space. 
However,  clouds,  as  well  as  dry  air,  partially  absorb  the  radiation  and  rc-cmit 
longwave  radiation  back  to  the  ocean  surface.  Thus  the  net  upward  energy  loss 
by  longwave  radiation  at  the  ocean  surface,  Rh,  is  corrected  for  the  downward 
radiation  by  the  clouds  and  the  air.  From  longwave  radiation  data,  Budyko 
(1978)  derived  a  scmi-cmprical  formula: 

Rh  =  a-\  bT -(a,  +/*,)/»  (17) 

The  dimensional  coefficients  a ,  /»,  ,  and  />,  arc 

a  =  -377.6  H'/ii"1,  l>  =  2.2  Wm~7K~\ 

=  -389.8  If7»»-2, />,  -  1.6 


7.  Basic  liquations  for  Perturbations 

When  the  coupled  system  is  perturbed  from  its  equilibrium  state,  the 
thermodynamic  feedback  mechanism  between  the  cumulus  clouds  and  the  oceanic 
mixed  layer  makes  the  perturbation  either  grow  (positive  feedback)  or  dampen 
(negative  feedback).  The  principal  purpose  here  is  to  study  thermodynamic 
feedback  mechanisms  between  clouds  and  the  oceanic  mixed  layer,  lienee  the 
energy  exchange  at  the  air-ocean  interface  is  a  primary  focal  point.  Therefore, 
wc  shall  neglect  Initially  the  perturbations  of  those  variables  which  arc  not  tli- 
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rcctly  related  to  the  exchange  at  the  air-ocean  interface.  From  the  basic  equations 
of  the  coupled  system  (R),  (10),  (12),  and  (13),  the  perturbations  satisfy  the  fol¬ 
lowing  equations: 


dri  \  dE  T  Mr  , 

dt  h  57*  w  i)n  ' 
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(Pr-E)S  „ 


at  hw 

8.  Three  Basic  Time  Scales 

Three  lime  scales  arc  found  from  (I8)-(2I): 

Cloud  time  scale 


hi 


h\. 


-I 


I  *Pr 


n  K  an 

Using  (9)  and  taking  hr  =  5cm,  we  have 

r„  ~  0.6  day 

SST  variation  time  scale 

I  BF 


-\  _ 
t7  = 


PSrnK  dT* 


(18) 


(19) 

(20) 

(21) 


(22) 


(23) 


(24) 


The  net  heat  loss  from  the  ocean  surface  is  computed  by  (14).  Assuming  that  0 
and  q  arc  determined  by  the  large-scale  atmospheric  motion,  we  have 

OF  ,»/,  .  CE  .  _ 

/V/»«”*(  b  _ t  r'  KC„ 
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CpnRjl  c" 


therefore 


i  P*cpavg  L2qt  CE 

Tv  =  1  - TT T-F-)ckch 


Pwfpw^w 


C/mKji  (  " 


(25) 


Cloud  and  SST  coupling  lime  scale 


(26) 


T-i  _  i  or  dE  .  tPr 

PwCpJ’w  <ln  dTw  <ln 
Substitution  of  (22)  into  (26)  lends  to 


Vr-[ - ,  ,  3"  2 

Pw^JlJcRJ* 


Using  (16)  and  (17)  wc  have 


(27) 


dF  *Wb  ~  Rs) 
dn  f)n 
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i 

Among  these  three  time  scales,  r„  is  (he  shortest.  The  other  two,  t,  and  t 
largely  depend  on  the  parameters  Cv  C,„  and  CK  ,  which  arc  functions  or  the  at¬ 
mospheric  stability.  Fig.4  shows  the  dependence  of  t7  and  r„  ,  on  the  atmospheric 
stability.  For  the  unstable  atmosphere  {hJL,  <  0),  t„7  ~  3-6  days,  and 
r7  ~20  -  30  days.  For  the  stable  atmosphere  (/>„//.„>()),  r„ —  0.3  —  I  years, 
and  t7  ~  I  —  3  years. 


Fig.4  Dependence  of  r  ,  and  t„  ,  on  the  atmospheric  stability 
parameter 


9.  Solutions 


Since  the  time  scale  for  cloud  feedback  is  so  much  shorter  than  that  for  SST 
feedback,  i.c.,  the  cloud  cover  perturbation,  /?',  almost  instantaneously  fol¬ 

lows  the  SST  for  the  temperature  feedback: 


MldTw 

dPr/dn  w 


(28) 


<V>,  , 

-r - 8.33  x  |()  1  ms  , 

(In 


Fq.(9)  shows  that 
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and  for  Tm  =  25 *C,  »/„.  =  0.lw/.«. 

CFua.-~~  ^  A.fti  x  IO_7/»/.v— 1  AT~ 

rhcrcforc,  a  I’C  change  in  SSF  implies  a  0.56  change  in  cloud  cover. 

Neglecting  the  small  terms  in  (he  prognostic  equations  for  7\,  and  S'  and 
eliminating  three  among  the  Tour  variables  /;',  S'  ,  and  from  (19),  (20), 
(21),  (2R),  we  obtain  the  second  order  differential  equation. 


(7V  2fi  -  y  ,  .  _K  d*  „Y  . 

**  -  —  Kr  +  *,  )-  +  K.r  +  tf  )  *  =  o  ,29) 

where  ^  represents  n  ,  7\,  .S', ,  and  h'^.  There  arc  two  nnndimcnsional  parame¬ 
ters, 

;|  Afjg(Pr  -  E)S 

+ *rl)  ’  V  ^V(vV+t7‘)  (30) 

that  indicate  the  relative  importance  of  mean  heal  and  mean  salinity  fluxes  in  the 
mean  surface  buoyancy  flux  (I  I).  The  general  solutions  of  (29)  have  the  form: 

tA  «  C,  cxp(<T| r)  I-  Cj  cxp(a2t)  (31) 

where  C,  and  C}  arc  integral  constants,  and  o,  and  a2  arc  the  eigenvalues  which 
are  the  rooLs  of  the  second-order  algebraic  equation 

0>  “  T=»r  (Vr  +  vV  +  -  "y  .  (vi  +  if1)2  -  0  (32) 

ir-p) 

10.  Instability  and  Oscillation  Criteria 

I  he  instability  criteria  for  the  thermodynamically  coupled  air-occan  system 


.-I  .  -K2 


<  0  decaying 

Rc(o)  •  =  0  neutral , 

>  0  growing  (33) 

where  a  is  or  the  root  of  the  second-order  equation  (32).  The  osc  ilia  tin  cri¬ 
terion  for  the  coupled  system  arc 


lm(a) 


=  0  nonoscil/atory 


i1  0  oscillatory 


The  roots  a,  and  a,  arc 


^2^  y 

"‘•2  =  2(y  -  /!)  ^  H  T7  ±  v/c  -  l'V2l{ 2/1  -  y)2 


(35) 
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The  condition  for  .he  general  nr  grnwmg/decnying  mndcs.  which  can  be  he- 
tlucctl  from  (38),  '»*'• 


2/i  -  y 

y-  H 


<0  decaying 
=  0  neutral , 


(36) 


>0  growing 

ami  .he  condition  for  oscillatory/nonnscillalory  mode*  is: 

>  n  nnnoscillatory 


(i  ~  )* 

r  | 

<p  oscillatory 

Separation  of  different  modern  ^’c  »  '  ^  o)"Xtoich  corresponds  to  the  western 
buoyant  damping  regime  (B  *  '  *’  J  .  .  _c  „f  plc  mean  surface  heal  and 

Pacific  warm  pool  regime,  the  re  altv controli„g  the  modes 
salinity  fluxes  in  the  mean  buoyarmy  flux  B  . « i  ^  mcMUrcd  by  y, 

of  the  coupled  air-ocean  system.  ...ji  factor),  however,  the  mean  up- 

St  the  same  results: 


V  > 


2/i 


. Nnnoscillatory  Damping 


2 /' 

<N* 

A 

p- 

A 

Oscillatory  Damping 

1  - 

2/i 

=*. Oscillatory  Growing 

2p>y> 

1+77 

2/i 

>y  >/‘ 

^Nnnoscillatory  Growing 

(38) 


1  + 


77 
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l;)f  '  Separation  of  different  model  in  the  (f,  Ml  plane 


1 1.  Conclusion 

The  feedback  between  the  cloud  and  OPBL  in  the  coupled  MABL  and  stable 
OPBL  system  is  investigated  by  a  coupled  model.  The  lime  scales  largely  depend 
on  the  stability  of  the  MABL.  For  the  stable  atmosphere,  the  two  time  scales  arc 
quite  long:  xT  ~  100-300  days,  and  r„r~  1-3  yr.  For  the  unstable  atm«»sphcrc, 
however,  the  two  time  scales  arc  much  shorter:  rr~  3-6  days,  anti  20-30 

days.  In  the  western  Pacific  warm  pool  regions,  the  MABL  is  usually  unstable. 
Therefore,  this  theory  may  provide  some  explanation  of  the  two  time  scales  (3-6 
and  20-30  days)  of  intense  convection  in  the  western  Pacific. 

In  the  coupled  system,  the  fresh  water  influx  at  the  ocean  surface  due  to  the 
excess  precipitation  over  evaporation  is  a  damping  factor  (negative  feedback). 
However,  the  surface  cooling  is  a  forcing  factor  (positive  feedback).  The  relative 
strengths  of  these  two  surface  fluxes  determine  the  mode  type:  decaying  or  grow¬ 
ing,  oscillatory  or  nonoscillatory. 
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